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Abstract 

The CP-violating effects to the subprocess gg — > H + H~ is studied in the mSUGRA 
scenario at the CERN LHC, by taking into account the experimental bounds of elec- 
tron and neutron EDM's. The CP violation effects in this process arc related to the 
complex phases of /i and Af in the mSUGRA scenario. In our calculation we consider 
small CP phases of /i and A f and neglect the effects of neutral Higgs bosons mixing. 
In this case the CP effects to the process mainly come from the complex couplings 
of Higgs-squark-squark. We find a strong dependence of charged Higgs boson pair 
production rate on the complex couplings in the MSSM parameter space. 

PACS: 14.80. Cp, 12.60.Jv,12.15.Ji 

Key words: Charged Higgs Boson, mSUGRA, CP Violation 



*This work was supported by National Natural Science Foundation of China. 



1 



I. Introduction 

One of the basic questions of particle physics is the nature of electroweak symmetry 
breaking. Various theoretical insights suggest the extensions from the minimal Higgs 
sector with one Higgs doublet, for example, the Minimal Supersymmetric Standard Model 
(MSSM). The MSSM predicts the existence of three neutral and two charged Higgs bosons 
h°, H°, A and ff^Q, in which two Higgs doublets must be introduced to provide masses 
to down- and up-type fermions through the superpotential and to keep the theory free of 
anomalies. The lightest neutral Higgs boson in the MSSM may be difficult to distinguish 
from the neutral Higgs boson of the standard model(SM). The exploration of additional 
Higgs bosons will be then very important to confirm extended Higgs sector. Discovery of 
a charged Higgs boson, H^, shows directly these extended versions of the Higgs sector. 
Beyond the discovery of Higgs bosons, reconstructing the Higgs potential will be necessary, 
and that requires the experimental study of the self-couplings of Higgs bosons. 

At the Large Hadron Collider(LHC), searching for these extra Higgs bosons is also one 
of the most important tasks. In much of the parameter space preferred by the MSSM, 
namely m H ± > m\y and 1 < tan j3 < mt / mi, , the LHC will provide several mechanisms 
to search for charged Higgs boson. A copious source of charged Higgs bosons with fair light 
mass are decay of top quarks t — > bH + [[| , since top quarks are produced with very large 
rates at the LHC. For a heavier charged Higgs boson the dominate production process 
is gb — > tH~&. The pair charged Higgs bosons can be produced via quark-antiquark 
annihilation [g] and gluon-gluon fusion^, [?], ||, [|]. Of course, Higgs pair production has a 
lower rate compared with the single production mechanism, but Higgs pair production via 
gg and qq annihilation involves the trilinear and quartic couplings of the physical Higgs 
bosons. So the measurement of these couplings will be useful in reconstructing the Higgs 
potential. 

The full MSSM Lagrangian consists of various soft masses and the higgsino /i parameter 



and has 45 CP-violating complex phases [10]. In principle, the CP phases may come 



from the soft SUSY-breaking gaugino masses M a (a = U(l)y , SU(2)l, SU(3)c), the soft 
bilinear Higgs mass mixing masses mf 2 ( or B^), the soft trilinear Yukawa couplings Af 
and the higgsino mass parameter ju. Assuming universality condition of the soft gaugino 
masses at the GUT (or Planck) scale, the gaugino masses have a common phase and 
the different soft trilinear Yukawa coupling are all equal at this scale. Due to the two 
global U(l) symmetries of the MSSM Lagrangian, the U(1)q Peccei-Quinn symmetry and 
the U(1)r symmetry acting on the Grassman-valued coordinates, not all CP phases are 
physical and the physical CP phases could be reduced to those of trilinear coupling Af 
and the fi parameter jjlj|], <J?^ and &A f , 

/i = He^, A f = \A f \j* A f. (1.1) 

The CP phases are constrained indirectly by the measurements of the fermionic electric 
dipole moment(EDM), mainly of electron EDM d e and neutron EDM d n . In order to 
satisfy the current experimental bounds of electron[12] and neutron [O] EDMs, there are 
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several possibilities to solve it. The first possibility is to let the CP phases very small, 
i.e.|$| < 10~ 2 [14]. The second possibility is that the CP phases are not very small 



but the first two generations of scalar fermions which contribute to EDMs are heavy [15]. 
The third possibility is that the partial cancellations among various contributions to the 
electron and neutron EDM's in a restricted parameter space and the CP phases do not 



have to be suppressed [|lq, 17]. 



In the Higgs sector of the MSSM, at the tree-level Higgs potential is CP invariant, since 
the mass parameters of the two Higgs fields and the vacuum expectation values tan (3 are 
real and \x enters only via |^| 2 . If one includes the radiative corrections involving soft- 
CP-violating Yukawa interactions in the theory, the tree-level CP invariance of the Higgs 



potential can be broken[18, It's so called explicit CP violation of the MSSM Higgs 
potential. Then the three neutral Higgs bosons can mix and none of the mass eigenstate 
Higgs boson is a CP eigenstate. For and \A$\ values larger than the arithmetic average 
of the scalar-top-quark masses squared, denoted as MgusY (^susy = ^-/^( m f 1 + m f 2 ))> 



the explicit CP-violating effects can be significant |18|[ . If and \Af\ values are smaller 
than Mg USY or the CP phases of \x and Af is much small, such CP-violating effects turn 
out to be negligible. 

Various SUSY processes including the effects of the CP phases have recently been 



considered, such as the effects in the neutral Higgs boson production at hadron collider [20], 
the effects in neutralino J2l| and chargino |22| production at LEP, at CERN LHC|^3), as 
well as at future electron-positron linear colliders pi]]. In this work we shall investigate 
the CP violation effects in process gg — > H + H~ in the context of the mSUGRA. In the 
numerical calculation we choose and \Aj\ smaller than Mg USY and the CP phases of 
\i and Af at the order O(10~ 2 — 10 -1 ) to neglect the effects of the neutral Higgs bosons 
mixing. While the constrains to the parameters from the measurements of the electron 
and neutron EDM's are taken into account. We organized the paper as follows. In Section 
II, the analytical calculations are presented. In section III there are numerical results and 
discussions and in section IV, a short summary is given. 

II. Calculation 

The Feynman diagrams contributing to the subprocess gg —* H + H~ and the total 
cross section of this subprocess can be found in ||. The effective Lagrangian for electric 
dipole interaction of a fermion can be written as 

C = -'-dfiio^F^. (2.1) 

At the one-loop level, the contribution in the framework of MSSM to the electron EDM 
comes from the virtual chargino and neutralino exchanges, for the neutron EDM it is 
contributed by the virtual chargino, neutralino and gluino exchange. By considering the 
electric dipole operator, the chromoelectric dipole operator and the purely gluonic six- 



dimension operator, one found [16] that there exist significant regions in the parameter 



space where cancellations between the gluino and the chargino exchanges reduce the elec- 



3 



trie and the chromoelectric contributions, and the purely gluonic parts lead to a dramatic 
lowering of the neutron EDM sometimes below the electron EDM value. In the case of the 
neutron EDM, the dominant chargino and gluino contributions appear with opposite sign 
over a large portion of the MSSM parameter space. In general, internal cancellations are 
more likely among the SUSY contributions to the neutron EDM than they are in the case 



of the electron. A recent calculation of Chang et al.[25] for the Barr-Zee-type[26] two-loop 
contributions to the EDMs, shows that the sbottom contribution is enhanced for a large 
tan /? f27fl . Therefore, in our analysis we choose suitable parameters and small CP phases 
to suppress the Barr-Zee-type twodoop contributions. The CKM CP phase in the SM can 
in principle contribute to the fermionic EDM, such as the electron EDM, but it turns out 
to be effective at threedoop level so that the contribution is very small. 

The mSUGRA is based on N = 1 supergravity grand unified theory, in which super- 
symmetry can be broken spontaneously via gravitational interactions and the electroweak 



symmetry is broken via radiative effect [28|. At the GUT scale, the CP-violating mSUGRA 
can be completely parameterized by just six input parameters, namely, nii/ 2 , mo, Aq, 
tan/3, and &a > where mi/ 2 , tuq an d Aq are the universal gaugino mass, scalar mass 
and the trilinear soft breaking parameter, $^ and are t ne CP phases of \i and the 
trilinear soft breaking parameter. From these six parameters, all the masses and coupling 
of the model are determined by the evolution from the GUT scale to the low electroweak 
scale. 

With the cross section of the subprocess, we can easily obtain the total cross section 
at hadron collider by folding the cross section of gg — > H + H~ with the gluon luminosity 



rl AC 

a{pp^ gg + X ^ H + H- +X)= d,T—^-a(gg^H + H-, at s = ts). (2.2) 

J4m 2 , Is dr 



I Am? . Is 

H+' 

where y/s and v§ denote the proton-proton and gluon-gluon c.m.s. energies respectively 
and — p£ is the gluon luminosity, which is defined as 

dC gg _ f 1 dx-i 



dr J T x\ 



fg{xi,<r)fg(—,Cf) 
•^1 



(2.3) 



Here we used r = x\x 2 , the definitions of x\ and x 2 are adopted from Ref.[^9[. In our 
numerical calculation we take the MRST(mode 2) parton distribution function f g {xi,Q 2 ) 
fl30| , and ignore the supersymmetric QCD corrections to the parton distribution functions 
for simplicity. The factorization scale Q was chosen as the average of the final particle 
masses mjj+. The numerical calculation is carried out for the LHC at the energy lATeV . 



III. Numerical Results and Discussions 

In this section, we present some numerical results of the CP effects to total cross section 
from the full one- loop diagrams involving virtual (s) quarks for the process pp — > gg + X — > 



4 



H + H~ + X. The input parameters are chosen as: m t = 173.8 GeV , mz = 91.187 GeV, 
nib = 4.5 GeV, sin 2 8\y = 0.2315, and a = 1/128. We adopt a simple one-loop formula for 
the running strong coupling constant a s 

aM = t^W^^W' (3 - 1) 

where a s (mz) = 0.117 and nj is the number of active flavors at energy scale \x. In our 
numerical calculation to get the low energy scenario from the mSUGRA, the renormaliza- 
tion group equations(RGE's) pl| are run from the weak scale Mz up to the GUT scale, 
taking all threshold into account. We use two loop RGE's only for the gauge couplings and 
the one-loop RGE's for the other super symmetric parameters. The GUT scale boundary 
conditions are imposed and the RGE's are run back to Mz, again taking threshold into 
account. The REG's for the gauge and Yukawa couplings, and the diagonal elements of 
the sfermion masses and gaugino masses are real. The phase of \x does not run because it 
cancels out in the one loop renormalization group equation of \i, and the magnitude and 
the phase of Aj do evolve [lq]. The decay widths of the intermediate Higgs bosons are 



considered at the tree level and these formula can be found in ref.[32]. 

In order to describe the CP violating effects in the process, we introduce the notation 
5 defined as 

where ctq is the cross section of this process with CP conservation and acp is the cross 
section with CP violations. In Fig.l we present the S versus CP phases of \x and A§. 
The input parameters are chosen as mo = 100 GeV,m\/2 = 150 GeV, Aq = 300 GeV , 
tan /3 = 4 and the CP phases of fi and Aq run from zero to 0.1-7T. When the CP phases are 
set zero, the masses of Higgs bosons are m h o = 97 GeV, m H o = 258 GeV, m A o = 254 GeV 
and niff+ = 266 GeV and the cross section of this process is 1.2 fb. In our numerical 
calculation we let = when §a runs from zero to 0.17T and vice versa. In this case, 
the CP phase corrections to the masses of Higgs bosons and squarks are the order of just a 
few percent. So the CP effects to the process at this case mainly come from the phases in 
the Higgs-squark-squark couplings. The CP effect of ^ and $a to this process increase 
as the phase becoming larger. The 5 is positive for $^ running in this region and can 
reach 12% when the $^ = 0.1-7T, but for $a ? t ne 5 is negative and can reach —46% when 
the &a = 0.17T. From the figure we can see that the effect of $>a i s larger than that of ^ 
to this process. There is a strong dependence of the production rates of the pair charged 
Higgs on the complex couplings in such MSSM parameter space. 

In Fig. 2 we also present the 5 versus CP phases of fi and Aq. The input parameters 
are chosen as tan/3 = 30 and the others are the same as those in Fig.l. When the CP 
phases are set zero, the masses of Higgs bosons are m h o = 106 GeV, m H o = 175 GeV, 
m^o = 174 GeV and m H + = 194 GeV and the cross section of this process is 44 fb. 
We still let = when runs from zero to 0.1-7T and vice versa. In this case, the 
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CP phases corrections to the masses of Higgs bosons and squarks are very small and the 
CP effects to the process mainly come from the CP phases in the Higgs-squark-squark 
couplings. The CP effects of $ M and to this process increase as the phase becoming 
larger. Both the 5 of ^ and &a to this process are negative. For $^ the 5 can reach 
—5.7% when the $^ = O.Itt and for §>a can reach —12% when the <&a = O.ln. The 
dependence of the production rates on the complex couplings is also obvious. 



IV. Summary 

In this work we calculate the CP-violating effects to the process pp — * gg + X — * 
H + H~ + X in the mSUGRA scenario at the LHC collider, by taking into account the 
experimental bounds of electron and neutron EDMs. The CP violation effects in this 
process are related to the complex phases of fi and Af in the mSUGRA scenario. We 
choose small CP phases of these parameters to neglect the mixing of neutral Higgs bosons. 
In this case the CP effects to the process mainly come from the complex couplings of Higgs- 
squark-squark. We find a strong dependence of the production rates of the pair charged 
Higgs boson on the complex couplings in the MSSM parameter space we chosen. 
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Figure Captions 

Fig.l The S versus CP phases of /x and Aq in mSUGRA with the parameters uiq = 
100 GeV, my 2 = 150 GeV, A = 300 GeV, tan/3 = 4. The solid curve is for <£> Ao running 
from zero to 0.17T while the phase of \i is set zero and the dashed curve is for $^ running 
from zero to 0.1 tt while the phase of Aq is set zero. 

Fig. 2 The 5 versus CP phases of ^ and Aq in mSUGRA with the parameters too = 
100 GeV, mi/2 = 150 GeV, A = 300 GeV, tan/3 = 30. The solid curve is for <$> Ao running 
from zero to 0.17T while the phase of \i is set zero and the dashed curve is for $^ running 
from zero to 0.1 n while the phase of Aq is set zero. 
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Fig-(1) 
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Fig. (2) 
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